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Multi-phase  ignition  of  hypergolic  propellants  Mono-methyl  Hydrazine  (MMH)  and 
Red  Fuming  Nitric  Acid  (RFNA)  is  studied  numerically  to  understand  fundamental  pro¬ 
cesses  such  as  gas  phase  ignition,  vaporization  and  liquid  phase  chemistry  for  characterizing 
ignition.  Such  understanding  will  be  critical  for  future  design  efforts  targeting  rapidly  re¬ 
peatable  cyclic  ignition  of  these  propellants.  Three  test  cases  are  considered:  gas  and 
liquid  phase  autoignition,  an  opposed  jet  diffusion  flame  and  a  liquid  opposed  jet  diffusion 
flame.  For  the  first  case,  three  reduced  chemical  kinetics  mechanisms  are  used  to  study 
autoignition  behavior  in  the  gas  phase  and  with  premixed  liquids.  In  the  second  case,  a 
laminar  diffusion  flame  of  the  gas  phase  reactants  under  varying  strain  rates  is  studied. 
The  third  case  investigates  the  interface  behavior  for  liquid-liquid  contact  and  its  effect  on 
gas  phase  ignition. 
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Nomenclature 

Mono-Methyl  Hydrazine  (CH3NHNH2) 

Red  Fuming  Nitric  Acid  (88%HN03,  10%N2O4  and  2%H20) 
Oxidizer  to  Fuel  ratio  by  weight 
Time,  s 

Molecular  viscosity,  Pa  s 
Mole  fraction 
Mass  fraction 

Thermal  conductivity,  W /  m  K 
Molar  mass,  kg/  kmol 
Molecular  diffusivity,  m2  /  s 
Density,  kg/  m3 
Pressure,  Pa 
Temperature,  K 
Enthalpy,  J /  kg 
Entropy,  J /  kg  K 

Production  (destruction)  rate  of  a  species,  kmol/  s 

Reynolds  Averaged  Navier-Stokes  equations 

Detached  Eddy  Simulation 

Large  Eddy  Simulation 

General  Equation  and  Mesh  Solver 

Species  number 
Reaction  number 
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I.  Introduction 


The  ignition  behavior  of  hypergolic  propellants  is  of  importance  for  rocket  applications  that  require 
multiple  restarts  and/or  thrust  control.  Additionally  the  high  specific  impulse  and  low  energy  requirement 
for  storage  at  atmospheric  conditions  make  these  propellants  a  prime  candidate  for  small,  maneuverable 
rocket  applications.  One  such  set  of  propellants  is  monomethyl  hydrazine  (MMH)  and  red  fuming  nitric 
acid  (RFNA)  which  are  liquids  at  atmospheric  conditions.  The  ignition  delay  of  these  propellants  is  an 
important  parameter  and  its  analytical  understanding  is  desirable  for  developing  new  designs  as  well  as 
advanced  control  strategies.  Such  understanding  requires  multi-physics  modeling  which  in  turn  depends 
upon  the  knowledge  of  the  sub-phenomena  involved  in  the  ignition  process.  The  goal  of  present  work  is 
to  characterize  the  underlying  physics  and  to  develop  a  modeling  methodology  for  investigating  hypergolic 
ignition  of  MMH  and  RFNA. 

Early  experiments  with  the  hydrazine  family  of  propellants  provided  system  level,  empirical  understanding 
of  macroscopic  effects  [1-5].  Saad  et  al.  [6]  separated  the  processes  involved  as  shown  in  Figure  1  and 
studied  some  of  the  sub-phenomena  independently.  They  reported  that  the  competing  processes  of  liquid 
vaporization,  reaction  and  the  effect  of  pressure  on  each  are  important  for  hypergolic  ignition.  Recent 
studies  [7-9]  confirm  that  the  hypergolic  ignition  of  MMH  and  RFNA  involves  complex  sub-phenomena: 
spray  formation  upstream  of  the  impingement  point,  instabilities  in  the  sheets  resulting  in  intermittent 
liquid  contact  and  aerosol  formation,  gas  layer  development  between  the  two  liquid  sheets  due  to  liquid 
reaction  and  vaporization  as  well  as  satellite  drop  interactions  and  micro-explosions  in  the  outer  extremities 
of  the  liquid  sheets.  The  rate  and  extent  of  each  sub-phenomena  determine  the  ultimate  outcome  for  practical 
systems,  in  terms  of  one  discernible  quantity,  namely  the  ignition  delay  [10].  Figure  2  shows  a  schematic 
representation  of  some  of  these  sub-phenomena. 


Figure  1:  Various  processes  as  described  and  studied  by  Saad  et.  al.  [6] 


Historically,  the  hypergolicity  of  this  family  of  propellants  was  attributed  to  thermal  activation  [11] 
and  the  resulting  gas  intermediates  were  noted  to  be  important  [12].  High  speed  measurements  obtained 
recently  [7]  confirm  that  the  gas  intermediates  are  at  a  temperature  of  550  A.  The  compression  within 
the  generated  gas  layer  was  thought  to  influence  the  ignition  [13]  and  the  associated  delay.  Formation 
and  decomposition  of  solid  intermediates  was  also  reported  [14-16]  and  were  attributed  to  be  the  cause  of 
pressure  spikes  during  re-ignition.  The  gas  phase,  known  to  contribute  a  large  share  of  the  total  heat  released, 
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Figure  2:  Sub-phenomena  involved  in  the  ignition  of  impinging  jets  of  MMH  and  RFNA 


showed  multiple  flame  regions  [11,17]  and  was  addressed  by  various  experiments  [18-20]  to  develop  reaction 
mechanisms  [21,22].  Recent  studies  have  improved  our  knowledge  of  the  thermodynamic  properties,  reaction 
rates  [23-27]  and  reaction  pathways  that  are  of  relevance  to  the  hypergolic  ignition.  However,  experimental 
evidence  points  to  the  gas  phase  ignition  being  contingent  upon  the  liquid  interactions  [7, 10, 28]  and  hence 
the  overall  effect  must  be  considered  for  the  accurate  prediction  of  ignition  delay. 

The  time  and  length  scales  of  the  sub-phenomena  involved  can  pose  a  significant  challenge  for  experi¬ 
mental  measurements  [9,29]  and  high-fidelity  physics-based  simulations  can  significantly  benefit  the  under¬ 
standing.  However  the  use  of  full  chemical  kinetics  mechanisms  in  addition  to  the  fluid  dynamics  can  incur 
significant  computational  costs  due  to  the  large  number  of  species  and  due  to  the  stiffness  introduced  by 
the  time  scale  differences  in  the  chemical  reactions.  A  tangible  reduction  in  this  expense  can  be  achieved  by 
reduced  chemical  mechanisms.  Three  such  mechanisms  for  MMH- RFNA,  due  to  Labbe  et.  al.  [30,31]  are  a 
recent  advance  and  need  to  be  studied. 

As  such,  hypergolic  ignition  phenomenon  has  been  an  object  of  study  for  more  than  four  decades  and 
many  of  its  aspects  have  been  investigated  experimentally.  However  detailed  understanding  of  hypergolic 
ignition,  which  is  a  cumulative  effect,  is  limited  due  to  relatively  less  efforts  on  the  modeling  front.  In  order 
to  remedy  this  and  use  substantial  advances  in  computational  methodology,  there  is  a  need  to  establish  a 
modeling  framework  for  hypergolic  systems.  Following  in  the  footsteps  of  experimental  investigations,  gas 
phase  kinetics  should  be  the  first,  followed  by  gas-liquid  interactions.  The  overall  objective  of  the  paper 
is  thus,  to  investigate  isolated  sub-phenomena  as  well  as  their  coupled  impact  on  the  ignition  delay  while 
providing  a  characterization  of  the  three  reduced  mechanisms. 

The  present  article  starts  with  a  brief  summary  of  the  CFD  solver  and  the  chemical  mechanisms  used  in 
this  study.  We  then  consider  a  suite  of  test  problems,  each  of  which  focuses  on  a  specific  aspect  of  the  full 
model.  The  test  problems  include  (1)  autoignition  to  study  the  gas-phase  chemical  kinetics  and  the  combined 
effect  including  reacting,  vaporizing  liquids;  (2)  an  opposed  jet  diffusion  flame  for  understanding  the  kinetics 
and  transport  in  a  diffusion  flame,  and  finally  (3)  an  opposed  liquid  jet  case  targeting  the  impingement 
region  and  intervening  gas  layer.  Results  and  discussion  of  each  part  are  summarized  at  the  end  providing 
a  direction  for  future  efforts. 
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II.  Model 


The  solver  GEMS  (General  Equations  and  Mesh  Solver  with  multiple  approaches),  used  in  the  present 
work,  is  a  finite  volume,  time  accurate  solver  capable  of  handling  a  generalized  fluid  (caloric  and  state 
equations).  A  hybrid  LES-RANS  (Large  Eddy  Simulation  and  Reynolds  Averaged  Navier-Stokes  equations) 
approach  is  used  for  turbulence  modeling  with  Wilcox’s  k  —  uj  model  for  the  RANS  regions  and  a  DES 
(Detached  Eddy  Simulation)  approach  in  the  resolved  regions  of  the  flow.  It  has  an  approximate  Reimann 
solver  with  comprehensive  preconditioning  and  employs  the  line  Gauss  Seidel  algorithm  for  the  solution 
of  the  linear  system.  A  multi-block  unstructured  mesh  is  handled  in  a  parallel  fashion  for  two  as  well  as 
three  dimensional  solutions.  For  modeling  the  liquid  phase,  the  homogeneous  equilibrium  model,  based  on 
Amagat’s  law  is  used,  whereby  the  total  volume  of  a  fluid  is  equal  to  the  sum  of  partial  volumes  of  its 
components.  The  liquid  phase  is  then  included  as  a  fluid  species  with  a  density  equal  to  the  liquid  phase 
density.  This  approach  has  been  successfully  demonstrated  in  previous  work  for  modeling  of  cavitating 
flows  [32] .  The  reader  is  referred  to  prior  articles  for  further  details  of  the  solver  [32-36] . 


A.  Equations  of  Motion 

The  equations  of  motion  can  be  written  as: 


dQ  dE^  _  dFvi 

dt  dxi  dxi 


H  +  Q 


(i) 


where  the  vector  Q  is  a  set  of  conserved  variables,  E{  and  Fv{  are  the  inviscid  and  viscous  flux  vectors 
respectively  and  the  vector  H  contains  the  source  terms  other  than  in  the  species  equations.  The  last  vector, 
D,  contains  the  sources  (or  sinks)  for  the  species  and  is  calculated  using  a  kinetics  mechanism  where  the 
reaction  rates  are  given  by  the  frequently  used  Arrhenius  expression  [37], 


Kf  =  ATbe-^r 


(2) 


where  kf  is  the  forward  rate  of  reaction,  A,  b  are  constants,  T  is  temperature,  Ea  is  activation  energy  and  R  is 
the  universal  gas  constant.  Accurate  reaction  rate  determination  can  require  additional  details  such  as  third 
body  efficiencies  and  low  pressure  limit  rates  [38,39].  At  low  pressures,  two  sets  of  Arrhenius  parameters  are 
typically  specified  where  the  rate  ko  is  effective  at  low  pressures  and  k0 c  at  high  pressures.  For  the  transition 
to  the  high  pressure  limit  (“fall-off”  region),  the  reaction  rate  can  be  obtained  by, 

k = (ttk)  f  <3) 


where  Pr  is  the  reduced  pressure  scale,  given  by 


MM  (4) 

^OO 

and  F  is  the  broadening  factor,  obtained  either  from  the  Troe  parameters  [38]  or  the  Tsang  and  Herron  [40] 
linear  fit.  In  Troe  form,  F  is  given  as, 


logF 


log  (Pr)  +  c 

2" 

_n  —  d  (logPr  +  c)  _ 

log  (Fcent) 


where  Fcent  is, 

Fcent  =  (1  —  a)  +  ae^~T^t2^  +  e^~t3^T>} 

and  cqti,  £2  and  £3  are  Troe  parameters.  Other  constants  in  the  calculation  of  Fcent  are, 


(5) 

(6) 


c  =  -0.4  -  0.67  log  (Fcent) , 

n  =  0.75  -  1.27  log  (Fcent)  (7) 

d  =  0.14. 
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For  the  Tsang  and  Herron  fit,  F  is, 


F  =  tsai  +  tsci2  x  T 


(8) 


where  tsa\  and  tsa 2  are  rate  parameters.  Once  the  forward  rate  is  determined,  the  reverse  reaction  rate  can 
be  obtained  from  equilibrium  constant  Keq. 

Thermal  and  transport  properties  for  each  species  are  calculated  using  the  curve-fit  tabulations  due  to 
McBride  et  al.  [41]  with  the  exception  of  the  molecular  diffusivity,  which  is  calculated  using  Chapman- 
Enskog  theory  [42] .  Mixture  thermal  properties  are  calculated  as  mass  weighted  averages  while  the  mixture 
transport  properties  are  found  using  mixing  rules  due  to  Hirschfelder  and  Wilke,  modified  by  Bird  [43,44] 
and  Mathur  [45].  According  to  these  modifications,  the  molecular  diffusivity,  of  any  species  j  into 

the  mixture  is  given  as, 


D 


i-V 


i—mix  —  _ _ _  ]\r 


(9) 


Mixture  viscosity  is  calculated  using  the  relation: 


N 


Xj  fij 


j= 1  2^1=1  *19  ji 


where  the  function  is  defined  as, 


-*('*es)"8B")  (•*») 


-1/2 


The  thermal  conductivity  of  the  mixture  is, 


x~l  (Taa+a1  j 

\3  =  1  2^1= 1  A,  J 


(10) 


(11) 


(12) 


B.  Chemical  Kinetics  Mechanisms 

The  reaction  mechanism  due  to  Anderson  [22]  contains  81  species  and  513  reactions  and  can  incur  significant 
computational  cost  for  transient  reacting  flow  computations.  Two  reductions  of  this  mechanism,  RCheml 
and  RChem2  are  considered  in  the  present  work.  The  first  mechanism,  RCheml  consists  of  25  species  and  98 
reactions  while  the  second  reduction  RChem2  includes  29  species  and  120  reactions  [30,31].  MMH  oxidation 
chemistry  proceeds  through  successive  hydrogen  abstractions  to  form  HONO  and  daughter  species  CHXN2. 
The  development  of  these  mechanisms  is  recent  and  the  original  as  well  as  the  two  reduced  sets  RCheml 
and  RChem2  consider  hydrogen  abstraction  first  from  the  nitrogen  atom  attached  to  the  methyl  radical, 
followed  by  the  other  nitrogen  atom  which  is  then  followed  by  scission  of  the  methyl  radical.  A  third  reduced 
reaction  mechanism  RChem3  comprising  41  species  and  200  reactions  is  a  latest  set  of  reactions  for  these 
propellants  [46].  Based  on  the  morpholine  reaction  mechanism  due  to  Li  et.  al.  [26],  it  allows  for  different 
hydrogen  abstractions  from  MMH,  including  one  from  methyl  radical.  Additionally  it  contains  species  that 
are  observed  to  form  an  aerosol  after  the  initial  liquid  interactions,  thus  facilitating  simultaneous  modeling 
of  the  liquid  and  gas  phase  chemistry.  A  broad  summary  of  the  three  mechanisms  is  provided  in  table  1. 


Table  1:  Reduced  chemical  mechanisms  summary 


Mechanism 

Species 

Reactions 

RCheml 

25 

98 

RChem2 

29 

120 

RChem3 

41 

200 
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C.  Liquid  Reaction  and  Vaporization  Model 


The  reaction  of  liquid  MMH  and  RFNA  is  known  to  yield  an  aerosol  at  an  elevated  temperature  of  550  K 
and  higher.  The  major  species  measured  by  Wang  and  Thynell  [7]  and  by  Liu  et.  al.  [24]  are  CH30N02, 
N20,H2()  and  N2.  Based  on  these  products,  a  global  reaction  can  be  formulated  (equation  13)  with  a  heat 
of  combustion  of  8  MJ /kg  —  MMH  and  an  adiabatic  flame  temperature  of  1209  K. 


2  CH3NHNH2  +  4 HN03  — »  2  CH30N02  +  N20  +  5  H20  +  2  N2.  (13a) 

Recent  observations  [47,  48]  of  gas  evolution  from  a  drop  of  MMH  impinging  on  an  RFNA  pool  indicate 
a  delay  of  the  order  of  100  //s,  suggesting  that  the  liquid  reaction  is  complete  before  then.  Due  to  such 
short  time  scales,  this  interaction  has  been  challenging  to  probe  experimentally.  It  has  been  estimated 
that  ionic  reactions  may  precede  the  decomposition  of  the  intermediates  to  form  the  aerosol  and  hence  the 
activation  energy  for  the  global  reaction  is  expected  to  be  negligible.  Since  the  intermediates  are  known  to 
decompose  at  a  temperature  [9]  comparable  to  the  boiling  point  of  either  liquid,  the  temperature  dependence 
of  the  reaction  is  anticipated  to  be  insignificant,  therefore  allowing  the  use  of  collision  frequency,  A ,  as  a 
single  parameter  for  determining  the  reaction  rate.  Four  values  of  the  collision  frequency  and  corresponding 
temperature  transient  are  shown  in  Figure  3(a),  indicating  that  the  collision  frequency  of  1016  corresponds 
to  a  delay  consistent  with  the  observed  gas  evolution.  The  rate  of  the  liquid  reaction  is  therefore  given  as, 

kf  =  ATbe~Ea/RT  (14) 

A  =  1016;  b  =  0.0;  and  Ea  =  0.0 


Vaporization  is  modeled  using  the  Langmuir  expression: 


mv 


Spvref  2.303 R  (  rA  t  ) 

V2ttMWRT 


(15) 


where  S  is  the  surface  area,  pvref  is  reference  vapor  pressure  at  Trej,  T  is  the  current  liquid  temperature, 
MW  is  molecular  mass,  R  is  universal  gas  constant  and  A Hv  is  the  heat  of  vaporization.  The  surface  area 
of  the  liquid  in  the  homogeneous  equilibrium  model  relates  to  the  gradient  of  the  volume  fraction  of  the 
liquid,  which  in  turn  varies  as  the  inverse  of  the  cell  size,  1/Ax.  It  is  apparent  that  knowing  the  properties 
of  fluids,  this  expression  can  be  easily  cast  in  Arrhenius  form,  allowing  the  phase  change  to  be  modeled 
as  a  chemical  transformation  where  the  difference  in  the  enthalpy  of  the  phases  provides  the  expected 
endothermic  behavior.  Combining  the  vaporization  model  with  the  liquid  reaction  for  different  surface  area 
S  (Figure  3(b))  shows  that  the  time  scales  involved  are  substantially  different  and  for  a  premixed  liquid,  the 
vaporization  is  a  slow  process  compared  to  the  liquid  reaction. 
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Figure  3:  Transient  temperature  for  liquid  reaction  and  vaporization,  O/F  =  1.5 


III.  Results  and  Discussion 

In  order  to  develop  a  systematic  framework  for  a  time  accurate  reacting  flow  simulation,  we  begin  with 
autoignition  of  the  premixed  reactants.  Due  to  their  importance  in  the  vaporizing  spray  environment, 
transport  phenomena  are  investigated  next  using  the  opposed  diffusion  flames.  Lastly  the  case  of  opposed 
liquid  jets  of  MMH  and  RFNA  is  used  to  explore  combined  effects  of  liquid  and  gas  chemistry  along  with  fluid 
dynamics.  This  set  of  cases  provides  a  foundation  for  subsequent  modeling  efforts  targeting  more  realistic 
experimental  configurations. 

A.  Autoignition 

The  constant  volume  or  constant  pressure  “bomb”  reactor  can  be  modeled  as  zero  dimensional  cases  in 
order  to  explore  the  transient  behavior  of  premixed  propellants.  The  initial  condition  is  set  with  a  given 
temperature  and  composition  at  a  pressure  of  one  atmosphere  and  the  transient  solution  is  sought  using  the 
commercial  solver  Cosilab.  In  addition  to  understanding  the  transient  behavior,  performance  of  the  reaction 
mechanisms  is  also  of  importance.  Particularly  relevant  for  a  multi-dimensional  time  accurate  simulation  are 
the  initial  temperatures,  mixture  ignition  limits  and  time  scales.  Validation  of  a  zero  dimensional  version  of 
the  GEMS  solver  is  also  carried  out  by  comparing  against  the  transient  results  from  Cosilab. 

The  main  parameters  of  the  reactor  problem  are  the  oxidizer  to  fuel  ratio  (O/F  ratio)  and  the  initial 
mixture  temperature.  Wang  et.  al.  [7]  found  that  the  aerosol  formed  due  to  liquid  interactions  can  reach 
temperatures  of  550  K  and  ignition  in  the  gas  phase  is  observed  beyond  this  temperature.  Figure  4  shows 
representative  results  for  both  constant  pressure  and  constant  volume  bomb  cases.  As  seen  from  Figure  4(a), 
a  mixture  temperature  in  excess  of  650  K  leads  to  ignition  delays  consistent  with  experimental  observations 
and  is  close  to  the  reported  value  [49].  Ignition  delay  predictions  by  RCheml  and  RChem2,  shown  in 
Figure  4(a),  are  close  while  third  mechanism  RChem3  predicts  a  smaller  ignition  delay  at  a  lower  mixture 
temperatures  and  has  a  lesser  slope  than  RCheml  and  RChem2.  For  a  stoichiometric  mixture,  predicted 
ignition  delays  match  at  700  K  for  RChem2  and  RChem3  and  at  730  K  for  RCheml  and  RChem3.  In 
comparison  with  constant  pressure  case,  transient  mixture  temperature  and  pressure  for  constant  volume 
case  are  higher  for  the  same  heat  release  and  therefore  a  lower  ignition  delay  is  to  be  anticipated.  Comparison 
of  stoichiometric  mixtures  at  various  initial  temperatures  (Figure  4(b))  shows  that  the  constant  volume  case 
has  a  similar  trend  as  the  constant  pressure  case  with  a  factor  of  two  lower  ignition  delay  values. 

A  typical  temperature  transient  for  these  propellants  is  shown  in  Figure  5(a)  for  all  three  reduced  mech¬ 
anisms.  It  shows  evidence  of  three  steps  -  a  pre-ignition  first  step  followed  by  two  steps  during  ignition.  The 
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Figure  4:  Ignition  delay  predictions  for  stoichiometric  mixture 


pre- ignition  step  starts  with  a  sharp  rise  followed  by  a  gradual  increase  to  a  temperature  of  1000  A,  while 
the  second  step  takes  the  mixture  to  a  temperature  of  2500  K ,  and  the  last  step  completes  the  reactions  to 
near  equilibrium  condition.  As  was  reported  earlier  [50],  RCheml  predicts  a  lower  temperature  than  the 
equilibrium  and  the  third  transient  is  absent  for  RCheml,  delayed  for  RChem2,  while  in  case  of  RChem3, 
the  third  transient  is  coincident  with  the  second  transient.  Initial  MMH  oxidation  occurs  through  successive 
hydrogen  abstraction  by  N02  while  later  oxidation  occurs  through  the  action  of  hydroxyl  (OH)  radical.  As 
seen  from  Figure  5(b),  the  OH  mass  fraction  clearly  corresponds  to  the  third  transient,  which  is  absent  for 
RCheml,  delayed  for  RChem2  and  coincident  with  the  second  transient  for  RChem3.  This  behavior  is  akin 
in  part  to  the  hydrocarbon  chemistry  and  is  a  result  of  the  methyl  radical  in  the  MMH. 


(a)  Temperature  (b)  OH  mass  fraction 


Figure  5:  Transient  for  O/F  ratio  1.6 


In  order  to  understand  the  asymptotic  state  predicted  by  the  three  mechanisms,  species  compositions 
at  two  times,  20  [is  apart,  before  and  after  ignition  are  plotted  against  the  equilibrium  values  in  Figure 
6.  Figure  6(a)  shows  the  nitrogen  containing  species  and  Figure  6(b)  shows  the  hydroxyl  radical  and  the 
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products  of  complete  combustion.  Comparison  of  RCheml  and  RChem2  shows  that,  except  for  NO,  the 
nitrogen  containing  species  predicted  by  RCheml  move  towards  equilibrium  while  RChem2  predicts  that 
the  nitrogen  containing  species  are  near  equilibrium  composition  even  before  ignition.  Although  similar  to 
RCheml,  RChem3  differs  in  that  all  the  nitrogen  containing  species  move  closer  to  equilibrium  during  the 
ignition.  The  hydroxyl  radical  mole  fractions  predicted  by  RChem2  and  RChem3  are  noted  to  be  closer  to 
the  equilibrium  values  as  compared  to  RCheml  predictions  and  the  resulting  product  composition  is  also 
closer  to  equilibrium,  indicating  that  the  pathways  for  hydroxyl  radical  are  important. 


RCheml  RChem2  RChem3 


RCheml  RChem2  RChem3 


Figure  6:  Species  transient  for  O/F  ratio  1.6 


Exploring  RChem3  in  greater  depth  is  necessary  for  its  use  with  the  liquid  phase  to  isolate  the  effects 
of  liquid  reactions  and  vaporization  from  the  gas  phase  kinetics.  In  order  to  map  ignition  with  RChem3,  a 
range  of  O/F  ratios:  0.1  to  11.0,  corresponding  to  equivalence  ratios  25  to  0.25,  is  considered.  The  rich  end  of 
the  O/F  ratio  represents  a  fuel  droplet  vaporizing  in  an  oxidizer  rich  environment  whereas  for  lean  mixtures, 
a  vaporizing  oxidizer  droplet  is  enveloped  by  fuel  vapor.  Additionally,  the  initial  temperature  is  known  to 
depend  upon  the  liquid  reactions  and  is  varied  in  this  case  as  a  parameter.  Corresponding  autoignition  delay 
trends  are  shown  in  Figure  7.  A  clear  distinction  is  seen  from  fuel  rich  side  to  oxidizer  rich  side,  confirming 
historically  known  behavior  that  an  oxidizer  lead  has  a  lower  ignition  delay  Boorady:1967,lecourt2004mmh  . 
It  should  be  noted  that  the  equivalence  ratio  scale  is  not  linear  and  comparable  limits  on  the  O/F  ratio  scale 
are  from  one  to  eleven.  Within  these  limits,  the  fuel  rich  side  exhibits  two  orders  of  magnitude  difference 
in  the  ignition  delay  for  a  given  initial  temperature.  The  oxidizer  rich  side  shows  less  than  half  an  order  of 
magnitude  difference  for  the  same  temperature.  The  sharpest  increase  in  the  ignition  delay  is  three  orders  of 
magnitude  occurring  from  stoichiometric  condition  to  an  equivalence  ratio  of  0  =  3  at  500  K .  There  is  also  a 
difference  between  initial  mixture  temperatures  below  and  above  500  K.  Lower  temperatures  show  an  order 
of  magnitude  drop  in  ignition  delay  per  100  A  increase  in  initial  temperature  while  higher  temperatures  show 
less  than  half  an  order  of  magnitude  drop  for  the  same  increase  in  initial  temperature.  This  is  an  important 
characteristic  of  these  vapors  and  is  expected  to  influence  the  ignition  in  a  vaporizing  spray. 
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(a) 

Figure  7:  Variation  of  ignition  delay  with  O/F  ratio  and  initial  mixture  temperature 


B.  Premixed  Liquids 

Experimental  evidence  as  well  as  gas  phase  auto-ignition  behavior  suggests  that  liquid  reactions  provide  an 
initial  temperature  rise  for  vaporization  of  the  propellants  and  act  as  a  trigger  for  the  gas  phase  chemistry. 
Therefore  a  logical  next  step  is  to  consider  premixed  liquids  at  the  interface  of  the  two  propellants.  Temporal 
evolution  of  such  an  interface  will  be  rapid  and  highly  localized,  however  to  isolate  and  study  the  underlying 
transient  processes,  auto-ignition  of  the  liquid  propellants  in  a  homogeneous  system  can  be  useful.  In  this 
case,  premixed  liquids  with  a  given  composition  and  at  room  temperature  and  pressure  are  allowed  to  react 
and  vaporize  and  the  vapors  react  further  to  ignite  in  the  gas  phase.  Liquid  phase  global  reaction  is  used 
along  with  the  vaporization  model  and  gas  phase  kinetics  are  handled  using  RChem3,  which  includes  the 
products  of  liquid  reaction,  providing  a  path  to  final  products  of  the  combustion.  Parametric  variation  of 
the  initial  mixture  composition  is  used  for  mapping  the  ignition  limits  while  temporal  history  of  the  species 
and  temperature  illustrate  local  evolution  of  the  liquid  reaction  sites. 

The  transients  of  some  of  the  intermediate  species  are  shown  in  Figure  8.  Species  derived  by  successive 
hydrogen  abstraction  from  MMH  show  an  initial  decrease  at  a  temperature  of  800  K  and  the  radical  HNCO 
is  seen  to  increase  from  this  point  till  ignition.  Catoire  [9]  found  that  HNCO  can  be  one  of  the  major  species 
sustaining  MMH/RFNA  combustion.  Indeed,  it  is  seen  that  the  HNCO  peak  occurs  at  the  peak  temperature 
while  HONO,  which  is  known  to  occur  through  oxidizer  decomposition  and  successive  hydrogen  abstraction 
peaks  earlier.  Methyl  nitrate  CH30N02  decomposition  into  CHsO  and  N02,  known  to  be  endothermic,  is 
seen  before  the  first  drop  in  the  hydroxyl  (OH)  radical  which  takes  the  mixture  to  approximately  1000  K. 
Availability  of  MMH  vapor  and  its  successive  derivatives  depends  upon  vaporization  and  these  are  seen  to 
build  up  gradually,  slowing  the  ignition  from  less  than  0.2  ms  (Figure  7)  for  premixed  vapors  at  800  A  to 
1  ms  for  the  premixed  vaporizing  liquids. 

Figure  9  shows  temperature  and  MMH  vapor  transients  for  various  initial  compositions  constituting  a 
parametric  evaluation  of  the  O/F  ratio.  It  is  seen  that  the  ignition  limit  on  the  oxidizer  rich  side  is  lowered 
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Figure  8:  Transient  species  with  liquid  reaction  and  vaporization  for  O/F  =  1.5 


from  an  O/F  ratio  greater  than  11.0  to  9.0.  Ignition  with  vaporizing  liquids,  unlike  premixed  gases,  is 
not  an  instantaneous  phenomenon  but  a  cumulative  occurrence  of  pre-ignition  reactions.  This  leads  to  a 
continuous  temperature  increase  up  to  the  ignition  point  in  contrast  to  the  almost  constant  temperature  of 
the  pre- ignition  gaseous  mixtures.  At  the  end  of  ignitable  range  (r  =  0.15,  3.3)  temperature  transients  lack 
the  sharp  initial  rise  to  1000  K  and  take  two  orders  of  magnitude  longer  to  reach  the  ignition  point  where 
a  sharp  temperature  gradient  is  observed  in  contrast  to  a  smooth  ignition  for  near  stoichiometric  mixtures. 
Vaporizing  sprays  can  produce  this  entire  range  of  conditions  and  can  therefore  be  one  of  the  reasons  for 
sharp  pressure  increases  seen  in  practice. 

C.  Opposed  Diffusion  Flame 

Many  studies  of  hypergolicity  of  this  family  of  propellants  [7, 15,48,51]  have  focused  on  a  drop  test  where  a 
drop  of  the  fuel  falls  onto  a  pool  of  oxidizer  or  vice  versa.  A  common  observation  from  such  tests  is  the  gas 
evolution  at  the  interface  and  subsequent  gas  phase  ignition.  A  vaporizing  spray  at  the  extremities  of  the 
liquid  sheets  also  produces  rich  and  lean  pockets  of  gaseous  mixtures.  The  length  scale  of  such  pockets  can 
be  an  order  of  magnitude  larger  than  a  vaporizing  drop  and  thus  mixing  in  such  regions  will  be  governed 
by  diffusion.  Experimental  evidence  suggests  that  the  diffusion  flame  of  MMH  consists  of  two  distinct 
regions  [15,  25]  and  it  has  been  hypothesized  that  the  first  flame- front  is  due  to  MMH  decomposition  and 
the  second  flamefront  completes  oxidation.  Therefore,  a  key  step  to  understanding  the  physics  in  the  gas 
layer  separating  the  liquid  sheets  or  the  diffusive  mixtures  elsewhere,  is  the  study  of  an  opposed  jet  diffusion 
flame. 

A  Steady-opposed  diffusion  flame  at  different  strain  rates  is  considered  and  a  solution  over  a  domain 
width  of  10  cm  is  obtained  with  smoothness  insured  by  adaptive  gridding  based  on  a  gradient  parameter 
of  0.2  and  a  curvature  parameter  limited  to  0.5.  Consistent  with  previous  section,  we  use  a  fuel  and  the 
oxidizer  temperature  of  800  K  and  vary  the  inlet  velocity  for  the  desired  strain  rate.  Structure  of  the 
flame  is  observed  using  selected  species  along  with  the  normalized  reaction  heat  rate  (HRR)  for  two  of  the 
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(a)  Temperature  (b)  MMH  mass  fraction 

Figure  9:  Transient  behavior  of  premixed  reacting  and  vaporizing  liquids  as  a  function  of  the  0/F  ratio 


mechanisms,  RCheml  and  RChem3. 


Figure  10:  Normalized  heat  release  rate  and  species  for  two  strain  rates  predicted  with  RCheml 


The  flame  structure  predicted  with  RCheml  for  two  strain  rates  is  shown  in  Figure  10.  At  a  low  strain 
rate  of  40  s-1,  the  normalized  heat  release  rate  (HRR)  shows  four  distinct  peaks  in  the  flame  region.  The 
first  hydrogen  abstraction  from  MMH  is  known  to  be  endothermic  and  coincides  with  the  fuel  side  trough 
in  the  HRR  and  the  production  of  the  final  MMH  derivative  species,  CH3NN.  After  the  last  non-methyl 
hydrogen  abstraction  from  MMH,  the  decomposition  proceeds  through  scission  of  methyl  radical  and  this 
step  is  known  to  be  exothermic.  The  first  peak  HRR  observed  on  the  fuel  side  coincides  with  the  destruction 
of  CH3NN  and  corresponds  to  the  last  step  in  MMH  decomposition.  The  next  HRR  peak  is  coincident 
with  methyl  radical  destruction  which  is  followed  by  a  CO  peak  indicative  of  a  two  step  oxidation  process. 
The  last  peak  HRR  corresponds  to  H2  and  CO  oxidation  and  is  the  highest.  Although  HONO  formation  is 
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expected  to  be  seen  on  the  fuel  side  due  to  hydrogen  abstraction,  the  peak  mass  fraction  of  HONO  is  seen 
on  the  oxidizer  side  due  to  endothermic  dissociation  of  HNOs  to  OH  and  N02.  At  a  relatively  high  strain 
rate  of  1280  s_1,  the  four  parts  of  the  flame  observed  earlier,  identified  by  the  peak  heat  release  rate  merge 
to  form  two,  where  the  first  peak  in  the  HRR  is  coincident  with  the  CH3NN  peak  and  the  second  peak 
corresponds  to  combined  destruction  of  CH3NN  and  CH3  as  well  as  oxidation  of  hydrogen  and  CO.  Noting 
the  change  in  the  magnitude  of  the  normalized  peaks,  it  is  clear  that  the  major  part  of  the  heat  release  is 
due  to  the  second  HRR  peak. 

For  the  same  strain  rates,  predictions  with  RChem3  are  shown  in  Figure  11.  Lower  of  the  two  strain 
rates  shows  two  flame  zones  identified  by  peak  heat  release  rates,  consistent  with  the  previously  observed 
dual  flame  of  hydrazine  and  N204  [15]  and  MMH  RFNA  combustion  seen  during  experimental  studies  at 
Purdue  university.  The  fuel  side  peak  heat  release  rate  is  the  highest  and  clearly  corresponds  to  MMH 
decomposition.  Unlike  RCheml,  MMH  decomposition  through  successive  H  abstraction  takes  place  within 
this  region  and  CH3  mass  fraction  remains  nearly  constant  up  to  a  diminished  second  peak  where  it  gets 
destroyed  with  CO  mass  fraction  reaching  a  peak  value.  The  second  highest  heat  release  rate  corresponds  to 
oxidation  of  CO  and  H2  by  OH  and  N02.  In  this  case,  the  HONO  concentration  drops  significantly  in  the 
second  peak,  indicating  that  alternative  pathways  in  RChem3  for  MMH  decomposition  are  effective.  Also 
noticeable  are  the  two  endothermic  regions  on  the  oxidizer  side,  which  are  comparable  in  magnitude  to  the 
peak  HRR  and  can  potentially  alter  ignition  behavior  in  case  of  a  multi-dimensional  simulation.  At  a  higher 
strain  rate,  the  flame  zone  consists  of  a  single  prominent  peak  corresponding  to  oxidation  of  H2  and  CO  in 
addition  to  diminished  HRR  peaks  corresponding  to  MMH  decomposition,  CH3  destruction  and  conversion 
of  remaining  derivatives  of  MMH.  A  notable  difference  between  the  predictions  of  RCheml  and  RChem3 
for  both  strain  rates  is  the  absence  of  an  endothermic  region  adjoining  the  flame  on  the  fuel  side.  For  both 
strain  rates,  the  width  of  the  flame  predicted  by  RChem3  is  twice  that  of  RCheml,  suggesting  a  need  of 
uniform  higher  spatial  resolution  in  multi-dimensional  simulations. 


(a)  Strain  rate  =  40  s  1  (b)  Strain  rate  =  1280  s  1 


Figure  11:  Normalized  heat  release  rate  and  species  for  two  strain  rates  predicted  with  RChem3 


D.  Opposed  Liquid  Jets 

Liquid  sprays  and  injectors  have  been  previously  studied  [52-57]  but  there  has  been  limited  inclusion  of 
liquid-phase  reactions  in  modeling  efforts  targeting  hypergolic  ignition.  As  found  from  experimental  evidence 
and  theoretical  computations,  liquid  interaction  is  a  key  process  in  hypergolic  ignition  of  these  propellants. 
For  an  impinging  jet  injector,  there  are  two  regions  of  such  interaction,  the  first  being  the  impingement 
point  and  unstable  sheets  formed  thereafter  while  the  second  is  the  less  intuitive  spray  interaction  at  the 
extremities  of  the  liquid  sheets  where  photographic  evidence  [10]  suggests  that  droplet  collisions  lead  to 
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multitudes  of  micro-explosion  events  [28,58],  triggering  gas  phase  ignition.  While  the  later  interactions  are 
substantially  challenging  to  investigate,  prior  interaction  can  be  numerically  studied  for  liquid  contact,  extent 
of  liquid  reaction,  resulting  gas  layer  and  its  temperature.  It  is  well  known  that  liquid  diffusion  is  an  order 
of  magnitude  slower  than  gas  diffusion,  making  the  liquid  reactions  a  surface  phenomenon.  Mixing  through 
convective  transport  is  therefore  the  primary  avenue  of  liquid  interaction  which  in  turn  depends  upon  gas 
layer  formation  and  interface  dynamics.  Interdependence  of  heat  generation  and  gas  layer  formation  needs 
to  be  investigated  to  determine  conditions  in  the  gas  layer  separating  the  liquid  sheets  that  are  formed  upon 
impingement. 

The  impinging  jet  injector  is  one  of  the  favored  types  of  injectors  for  these  propellants  and  building 
towards  that,  two-dimensional  opposed  impinging  liquid  sheets  are  considered  first.  The  sheet  width  is 
assumed  to  be  0.25  mm  which  is  a  third  of  the  orifice  radius  of  the  injector.  The  span  separating  the 
impinging  jets  is  1  mm  and  open  to  atmospheric  condition.  Relatively  smaller  dimensions  are  specified  in 
this  case  for  managing  the  computational  expense  of  obtaining  a  grid  independent  solution.  The  behavior 
of  reacting  vaporizing  liquids  is  considered  without  the  additional  expense  of  gas  phase  chemistry.  Plug  flow 
of  liquids  at  300  K  is  specified  at  the  boundaries  so  as  to  have  equal  momentum  of  the  impinging  streams. 
Resulting  velocities  are  15  m/s  for  RFNA  and  20  m/s  for  MMH  due  to  the  density  difference  between  the 
two  fluids.  To  isolate  the  effect  of  the  gas-liquid  interface,  turbulence  is  not  considered.  For  grid  independent 
solution  of  the  gas  layer  between  the  liquids,  a  series  of  grids  is  chosen.  It  is  found  that  a  grid  of  cell  size  1  fim 
provides  sufficient  resolution  and  there  are  5  to  8  cells  within  the  gas  layer.  With  this  grid,  the  computation 
is  carried  out  and  transient  and  average  quantities  are  monitored. 
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Figure  12:  Resolved  transient  liquid-gas  interface 


A  series  of  transient  temperature  contours  (Figure  12)  shows  that  after  initial  planar  contact,  the  interface 
rapidly  becomes  corrugated,  resulting  in  rolling  up  of  liquids  at  multiple  sites.  Such  sites  react  violently  to 
produce  gas  bubbles  which  are  squeezed  out  by  the  incoming  liquid.  Although  initial  contact  and  mixing 
leads  to  temperatures  in  excess  of  500  A,  rapid  bubble  generation  prevents  further  contact  of  the  two  liquids 
and  the  gas  layer  cools  down  to  a  temperature  of  less  than  450  K  as  seen  from  averaged  results  (Figure  13) 
for  two  grid-sizes. 

Parametric  experimental  investigation  [59]  of  jet  diameter  and  velocity  have  shown  that  adequate  liquid 
contact  is  critical  for  repeatable  ignition  and  modeling  can  help  understanding  hypergolic  ignition  further, 
which  is  necessary  to  design  injectors  that  are  capable  of  rapidly  repeatable  cyclic  ignition.  It  should  be 
kept  in  mind  that  there  are  two  clear  distinctions  in  the  case  of  the  impinging  jet  setup  -  firstly,  three 
dimensionality  which  allows  for  additional  expansion  and  slowdown  of  the  sheets  and  secondly,  a  component 
of  momentum  along  the  jet  axis  which  causes  shear  layer  development  and  enhances  mixing.  In  the  region 
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close  to  the  impingement  point  however,  as  the  gas  layer  acts  to  support  transverse  momentum  of  the  jets, 
it  is  likely  that  rapid  contact  and  limited  heat  generation,  as  seen  from  this  case  will  prevail.  Future  studies 
will  focus  on  extending  these  results  to  three  dimensional  impinging  jet  configurations.  Preliminary  results 
for  gas-phase  propellants  are  presented  in  a  companion  article  in  this  conference.  Liquid-phase  computations 
will  be  undertaken  in  the  near  future  and  will  be  closely  tailored  to  current  experimental  studies  at  Purdue 
university. 


Figure  13:  Average  temperature  and  reactant  mass  fractions  at  the  liquid-gas  interface 


IV.  Summary 

Three  fundamental  problems  -  autoignition,  opposed  diffusion  flame  and  liquid  opposed  jets,  are  inves¬ 
tigated  for  understanding  the  ignition  phenomenon  of  hypergolic  propellants  MMH  and  RFNA.  Gas  phase 
auto-ignition  is  modeled  with  three  reduced  kinetics  mechanisms  and  shows  an  initial  temperature  require¬ 
ment  of  650  A.  With  respect  to  mixture  composition,  the  ignition  delay  is  observed  to  increase  sharply  for 
fuel  rich  mixtures  while  fuel  lean  mixtures  show  the  ignition  delay  to  be  relatively  independent  of  initial 
composition.  This  is  in  agreement  with  the  established  practice  of  providing  an  oxidizer  rich  environment  for 
ignition  of  MMH-RFNA  combinations.  Gas  phase  ignition  by  virtue  of  its  initial  temperature  requirement, 
depends  upon  liquid  reaction  for  activation  and  including  a  liquid  reaction  and  vaporization  shows  different 
behavior  than  the  premixed  gas  auto-ignition.  Unlike  the  premixed,  pre-heated  gaseous  reactants,  which 
showed  pre-ignition  reactions  at  an  almost  constant  temperature,  the  homogeneous  liquid  mixture  shows 
continuous  pre-ignition  temperature  rise,  recovering  the  step  like  temperature  transient  in  the  limit  of  rich 
and  lean  mixture  extremes.  In  case  of  a  practical  device,  operating  outside  the  continuous  temperature 
increase  regime  may  be  an  issue  for  achieving  repeatable  smooth  ignition. 

The  opposed  diffusion  flame  investigation  indicates  multiple  zones  of  heat  release  corresponding  to  various 
aspects  of  the  chemistry.  At  low  strain  the  zones  corresponding  to  MMH  decomposition  and  subsequent 
methyl,  CO  and  hydrogen  oxidation  are  distinct  and  have  comparable  heat  release.  High  strain  rate  on  the 
other  hand  shows  merging  of  several  regions,  particularly  on  the  oxidizer  side,  with  a  diminished  contribution 
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from  MMH  decomposition.  An  endothermic  region  exists  adjacent  to  the  oxidizer  side  and  can  act  as  a 
delaying  mechanism  for  the  ignition  of  the  diffusing  reactants. 

Opposed  liquid  jets  provide  a  further  insight  into  liquid  mixing,  reaction  and  dynamics  of  the  gas  layer 
trapped  between  liquid  sheets.  The  gas  layer  develops  as  initial  products  formed  from  liquid  reaction  cannot 
stably  support  the  momentum  of  the  liquid  jets,  leading  to  intermittent  contact  and  bubble  formation 
between  the  two  liquids.  Local  temperatures  can  exceed  the  lower  threshold  of  the  gas  phase  ignition  but 
average  temperatures  in  two  dimensional  sheets  are  seen  to  be  restricted  to  450  K.  Future  work  will  focus 
on  extending  these  results  to  practical  three  dimensional  experimental  configurations,  wherein  additional 
effects  may  contribute  to  the  occurrence  of  ignition. 
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